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ABSTRACT 

In this paper, we aim to estimate the vertical gradients in the rotational velocity of the Galaxy. 
This is carried out in the framework of a global thin disc model approximation. The predicted 
gradient values coincide with the observed vertical fall-off in the rotation curve of the Galaxy. 
The gradient is estimated based on a statistical analysis of trajectories of test bodies in the 
gravitational field of the disc and in an analytical way using a quasi-circular orbit approximation. 
The agreement of the results with the gradient measurements is remarkable in view of other more 
complicated, non-gravitational mechanisms used for explaining the observed gradient values. 
Finally, we find that models with a significant spheroidal component give worse vertical gradient 
estimates than the simple disc model. In view of these results, we can surmise that, apart from 
the central spherical bulge and Galactic halo, the gross mass distribution in the Galaxy forms a 
flattened rather than spheroidal figure. 
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1 INTRODUCTION 

Recently, apart from other local and global chara cteristics of the rota- 
tion speed of the Galaxy, such as rolling motion, ^ Levine et alj j2008l) 
determined the fall-off rate 7 = —22 ± 6 km kpc~^ in the rota- 
tion speed from the Galactic mid-plane. This gradient estimation was 
obtained by fitting a linear profilqj 



X 

%^Vy,{r,z) = v^{r,0) +71^1 +Sz 



(1) 



■ to rotation measurements in the vicinity of the Galactic mid-plane: 
l^l < 0.1 kpc, r G (3,8) kpc. 

The principal purpose of our work is to reconstruct the vertical 
gradient magnitude in a simple model of the Galaxy and to find out 
to what extent the gradient behaviour is dependent on the geometry 
of mass distribution. In particular, the gradient value can be very well 
reconstructed in the global thin disc model, provided the disc com- 
prises gross dynamical mass ascertained from the rotation curve of 
the Galaxy. 



1.1 Motivation of the present work 

There is a suggestive heuristics behind the above simple ansatz for 
the rotation speed. This naturally leads us directly to the global disc 



^ Throughout this paper we use cylindrical coordinates (r, ip, z) distin- 
guished by the galactic mid-plane 2 = and the axis of rotation r = 0. 



model approximation as a mean of determining the vertical gradient 
of the rotation speed in flattened galaxies, in particular, in the Galaxy. 

It is a simple matter to note that every axisymmetric function 
/(r, z) in the cylindrical coordinate frame can be written (for z ^ 0) 
as 

/(r, z) = f{r, 0) + r(r, z) \z\ + A(r, z)z (2) 
Here, F and A are z-symmetric functions defined as 

r(r, 2) = f(r,\z\) + f{r^-\z\)-2f{r,0) ^ ^(^^ ^-J ^ /('', N I ) - / (r, - I z I ) ^ 

Consider now the azimuthal velocity in place of the function /, 
/(r, z) = v^{r, z). We expect i'<^(r, 2) to be nearly 2-symmetric. In 
this case, the difference /(r, |2|) — /(r, — |2|) is small compared to 
the symmetric part f{r, \z\) + f(r, — \z\) ~ 2/(r, 2). From the La- 
grange mean value theorem, it then follows that F « = 
fxi^' 0\c=hzSgn{z), < h < 1 and A ~ 0. When, in addition, 
v^{r, 2) falls off almost linearly with distance off the mid-plane and 
does not change significantly with r in the considered region, then F is 
nearly constant. In this case, we can write (r, 2) ~ t;^(r, 0)+7 j2j, 
where 7 is a constant such that 7 ~ F. Thus, to the first order in 2, 
constant 7 measures the fall-off rate or the vertical gradient of the 
speed of rotation. In the same order of approximation, we can in- 
clude the correction from the antisymmetric part. Then (r, 2) ~ 
Vip{r,0) + 7 I2I + 5z, provided A is also almost constant in the con- 
sidered region. Under these conditions, we can use model ([TJ to es- 
timate the characteristic magnitudes of F and A. Constant 5 is then 
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called a rolling parameter, as \5z\ measures the difference of speeds 
on opposite sides of the galactic mid-plane. In this approximation, 
parameters 7 and 5 are global characteristics of the rotation speed in 
the considered region of coordinates r and z. Of course, the linear 
approximation can also be used when functions V and A are not al- 
most constant in this region. Then parameters 7 and 5 characterize, 
on average, the fall-off rate and rolling in this region. 

In the z-symmetric case, customarily assumed in simple mod- 
elling of mass distribution in galaxies, we necessarily have A = 0. 
Thus, the rolling parameter cannot be determined by such models, 
but the fall-off rate can be still determined. If we can show in such a 
galaxy model that F indeed is almost constant in the considered re- 
gion and its predicted value is correct, then we can regard the model 
of mass distribution as consistent to the first order of approximation 
with the method used to determine 7 from observations of rotation in 
the Galaxy based on the linear fit (T). The same applies when F is not 
almost constant; then, 7 characterizes, on average, the fall-off rate in 
the considered region. What counts is that the linear fit to observa- 
tions l[T) and a model value of the fall-off in a linear approximation 
are consistent with each other. 

In this paper, we attempt to realize the above programme of find- 
ing a simple model, which in the linear approximation would be able 
to predict the correct value of 7. Although expression for /(r, z) on 
the right-hand side of equation l|2j, which we used above as a start- 
ing point to show that approximation l[T) is well grounded, is always 
valid - has the same analytical properties as f{r, z)~ the approximate 
expression ([TJ is not analytical at 2; = as jzj cannot be differenti- 
ated. This is the same type of singularity that is encountered in the 
thin disc model in which the gravitational field is everywhere smooth 
apart from the plane z = 0\ the same goes for rotation on almost cir- 
cular orbits close to this plane. The thin disc model is thus naturally 
compatible with approximation l[T). Of course, mass distribution in a 
real flattened galaxy does not form an infinitely thin disc and {r, z) 
is a smooth function of which equation l[T) is a cusp-like approxima- 
tion. However, a more realistic model of mass distribution (say, finite 
width disc) in the same linear approximation should lead to similar 
results, although in an unnecessarily more complicated way. We have 
therefore decided to model vertical gradients in the rotational veloc- 
ity of the Galaxy in a global thin disc model approximation. Later, we 
compare predictions of the model with the observed gradient value. 
We also find from numerical simulations of test bodies in the gravi- 
tational field of the disc that the vertical gradient is almost constant. 
This again gives further support that our model is consistent with the 
gradient measurement method l[T) (constancy of the fall-off rate is 
also in accord with observations of the gradient in other galaxies). 

When the fall-off rate in the considered region is observed to 
be constant (which is the case for other galaxies), the important fea- 
ture of linear approximation is that, in order to model the fall-off 
rate 7 in this region, we can investigate the behaviour of the model 
speed not necessarily very close to the galactic mid-plane, but also at 
heights greater than the stellar disc width. In this situation, the scale- 
height should not be important as additional non-linear terms, which 
could be added to the fitting profile l[T), would be less significant. In 
other words, when the gradient value is observed to be constant, we 
should expect that models with intrinsic scaleheight (such as finite 
width discs) or without such scale (such as a thin disc) would predict 
com parable estim a tions f or the gradient value in the region. 

iLevine et a found the observed gradient values to be 

roughly four times greater than those predicted by a finite width expo- 
nential disc of stellar mass. They correctly concluded that to explain 
the gradient value some other process besides simple gravitational 
physics must have been taken into account. It is important to note that 
this situation cannot be remedied by taking into account other mass 
components, such as the dark halo, because close to the Galactic mid- 
plane any spheroidal mass component gives negligible contribution 
to the total gradient value. It vanishes on the mid-plane and only be- 



comes significant further away (this issue is illustrated in more detail 
in the next paragraph). This is the principal difference between contri- 
butions of a spheroidal and a flattened mass distribution to the overall 
gradient value close to the Galactic mid-plane. 

Interestingly, we show that the high gradient values could be ex- 
plained by simple gravitational physics if only the premise about mass 
distribution in Galaxy were changed. To this end, we have assumed 
that perhaps except for the central core, the whole dynamical mass 
forms a flattened, disc-like, rather than spheroidal, structure. This hy- 
pothesis enables us to use a global thin disc model for the approximate 
description of a flattened mass distribution comprising all dynamical 
mass (i.e. that inferred from Galaxy rotation). We can therefore say 
that high gradient values may indicate simply that gross dynamical 
mass forms a flattened mass distribution, and not that other processes 
besides simple gravitational physics would be needed to explain the 
gradients. 



1.2 Efficiency of flattened and splieroidal mass distribution in 
explaining rotation speed fall-off 

For illustration, let us compare vertical gradients for a disc-like and a 
spherically symmetric mass distribution with the same rotation law in 
the z = pla ne. For simplicity, we can consider the rotation law of a 
Kuzmin disc jSinnev & Tremain^l 19871) 
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It is a simple matter to show in the quasi-circular orbit approximation 
(we discuss this approximation in Section|2]l that the vertical gradient 
of azimuthal velocity is 
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However, the absolute gradient value would be much lower for a 
spherically symmetric model with the same rotation curve 
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and it would vanish at the disc plane. The gradients would become 
comparable only at some height off the disc. 

In general, the contribution to the vertical gradient of rotation 
from a spherical potential also vanishes in the galactic mid-plane, 
whereas a similar contribution from the gravitational field of a disc 
remains non-zero. From this, it also follows that to attain the actual 
gradient values at larger heights off the galactic mid-plane (outside 
the main concentration of masses, where the 'thin disc gravitation' is 
comparable with the 'wide disc gravitation'), the disc model gradient 
must decrease in absolute value, while the spherical model gradient 
must increase in absolute value. Qualitatively, this explains why the 
vertical fall-off of rotation is, on average, weaker when the gravity 
of a flattened, disc-like component is dominated by the gravity of a 
spheroidal component. The presence of a massive spheroidal com- 
ponent reduces the overall vertical gradient in two ways: (i) because 
(as we have already seen) the contributions to vertical gradients from 
spheroidal components are small; (ii) because the presence of a spher- 
ical component uses up masses that would otherwise be present in the 
disc component, enhancing the disc-like contributions to the overall 
gradient. 

1.3 Outline of the results 

In the framework of a global thin disc model approximation for the 
Milky Way, we use various methods for estimating the vertical gradi- 
ent in rotation speed. All the methods lead to comparable results. In 
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particular, we find the gradients for 0.22 ^ \z\ ^ 2.62 kpc by us- 
ing an averaging method, mimicking realistic gradient measurements 
in galaxies other than the Milky Way, and for < l^;] < 3.6 kpc 
by analysing the motion of test bodies in the gravitational potential 
of the Galaxy described in the thin disc model approximation. The 
gradient values are compared with the observed value determined by 
iLevine et al.l ( 1200 8l) from measurements of the rotation of the Milky 
Way in the region \ z\ < 0.1 kpc. Although our various gradient esti- 
mates are carried out at larger z, we consider this comparison possi- 
ble. One argument for this is that the value of the velocity fall-off is 
almost constant with the distance off the mid-plane we encountered in 
our analysis. Surely, it would be best to measure the gradient values 
in the Galaxy at heights larger than 0.1 kpc; however, such data are 
not available. However, we can utilize the observational fact that the 
gradient is constant in other galaxies. In addition, we also calculated 
the gradient analytically in the region \z\ < 0.1 kpc, again obtain- 
ing results consistent with the observed gradient value in this region. 
The gradient values we predicted us ing our simpl e model agree amaz- 
ingly well with those obtained by iLevine et al.l |200^ even though 
the Galaxy rotation curve we utilized is unrelated to that paper. This 
agreement is striking in view of the fact that the gradients are usually 
considered to originate from more complex physics than gravitational 
alone, to mention only radial pressure gradients, magnetic tension, or 
galactic fountains. 

For completeness, we compare our results with those predicted 
for the gradient by a customary three-component model or a maxi- 
mum halo model. It turns out that the gradient features are crucially 
dependent on whether the gross mass distribution is more flattened or 
more spheroidal, and we have already encountered this feature when 
considering a simple example in Section [L2l We find that the gradi- 
ent estimates for the Galaxy predicted in the disc model are in better 
agreement with observations than analogous predictions of models 
with a significant spheroidal mass component. 

Finally, we apply the disc model and estimate the vertical gradi- 
ents in the galaxy NGC 891, where rotation was measured far from the 
galactic mid-plane, again obtaining results consistent with measure- 
ments. These results show that the simple thin disc model approxima- 
tion performs well in reconstructing the vertical gradient properties. 
Of course, this does not mean that galactic discs are thin. However, it 
does provide a strong argument for the fact that gross mass distribu- 
tion in the Milky Way and in NGC 891 might be flattened rather than 
spheroidal. 



1.4 A model (SHO) rotation curve of Milky Way Galaxy 

Apart from the main reference jLevine et alj200^ . the observational 
basis for our analysis is a unified set of data representing the rota- 
tion velocity of the Galaxy collected bv ISofue. Honma & Omodakal 
(2009, hereafter SHO). Based on these data, we derive a substitute 
thin disc surface mass density accounting for the Galaxy rotation and 
representing a sort of projection of Galaxy masses on to the Galac- 
tic mid-plane. If gross Galaxy matter forms a flattened structure, as 
seen in other spiral galaxy pictures (apart from the galactic luminous 
halo and the central spheroidal bulge, the latter, nevertheless, com- 
monly described by an equivalent de Vaucouleurs disc), then the ac- 
tual Galaxy gravitational potential should be better approximated by 
the disc model than by a model with a significant spheroidal compo- 
nent comprising most of the dynamical Galaxy mass inferred from the 
Galaxy rotation. The data set is a compilation of several independent 
measurements of the rotation, suitably adjusted to each other. Unfor- 
tunately, the rotation data are poor at larger distances where the mea- 
surement points are scattered randomly and are determined with large 
errors. This leaves too much room for arbitrariness, as neither declin- 
ing nor even rising rotation could be excluded at outer radii (however, 
this inaccurately measured, remote Galaxy region is not very impor- 



tant to our analysis). Therefore, instead of the unified measurements, 
we decided to use a published model rotation curve that accounts for 
the data satisfactorily, mainly in the 3 — 8 kpc region in which we are 
interested. This curve was obtained by SHO for the unified measure- 
ment data by fitting a three-component model, including in addition 
a wavy ring pattern superposed on the exponential disc component. 
For brevity, we refer to this particular fit as the SHO rotation curve 
and denote it by Va ■ The SHO rotation curve agrees, as a constrained 
least-squares fit, with the Galaxy rotation curve almost perfectly out 
to large radii (however, any other fit, such as a linear combination 
of Hermite polynomials, a Fourier series, etc., anything one wishes, 
would also be acceptable, if only the fit was sufficiently smooth and 
close to the data in the sense of some norm, such as the least-squares 
method norm). The rings allowed the authors to reproduce the ob- 
served dips in rotation. Leaving the SHO mass model aside, the SHO 
fit alone can be treated as an actually measured rotation curve of the 
Galaxy, as it agrees very well with real rotation measurements in the 
internal Galaxy region, accounting for qualitative characteristics of 
the rotation. This curve agrees also with additional constraints im- 
posed on the curve in SHO and known from separate precise single 
rotation measurements[f| 



2 VERTICAL VELOCITY GRADIENTS OF AZIMUTHAL 
VELOCITY IN THE GLOBAL THIN DISC 
MODEL:QUASI-CIRCULAR ORBIT APPROXIMATION 

The surface mass density (T(r), representing a flattened mass distri- 
bution projected on to the Galactic mid-plane and corresponding to 
SHO rotation curve, is found in a global thin disc model approxima- 
tion from 



E ^] 



+ 



E ( — ^ 

X — r-' 



(3) 



We derived this equation in I JalochZ Bratek & Kutscherd §00i). In 
contrast to the equiva lent relation known from the classic textbook 
on galactic dynamics teinnev & Tremain3ll987h . it does not contain 
derivatives of rotation. Because no additional constraints on mass dis- 
tribution are taken into account in this paper, the integration in in- 
tegral (O must be cut off at some radius Rcr (when integration be- 
comes cut off, the resulting cut-off errors, as a result of some bound- 
ary effects, are different both in value and behaviour from those re- 
sulting in the classic formula). Ac cording to a criterion derived in 
iBratek. Jalocha & Kutscheral j2008l) , the uncertainty in (j{r) determi- 
nation, resulting from the cut-off, can be neglected when < | Ra . 
Here, the cut-off radius equals — 20 kpc, which is the radial ex- 
tent of Vcr{r). In particular, the criterion is satisfied inside the ring 
r G (3, 8) kpc, which is the region that interests US0 

Given a surface mass density, azimuthal velocities of test parti- 
cles in the gravitational field of a thin disc can be estimated based on 
the radial force equilibrium condition 



^ For clarity, it .should be noted that a square root in the rotation law for the 
assumed dark halo mass distribution was unwittingly missed in SHO. This 
is probably not a mere misprint as a plot of the dai'k halo contribution to the 
overall rotation curve presented in their paper is parabolic at the centre, despite 
being linear for such a halo. Obviously, this mistake is not important to our 
analysis and does not affect our results. 

The lower bound for the Galaxy mass in the global disc model, integrated 
out to cut-off radius Rcr, is 1.12 X IO^^Mq for the SHO rotation curve, using 
a substitute surface density defined in equation (5). 
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-gr{r,z), 



(4) 



Here, gr{r, z) is the radial component of tlie gravitational acceleration 
(this approximation is justified later), and hence 



= G 



(r - X cos (f) a (x) xdxdy 
(r^ + X^ — 2rx cos cp + z^'f'^'^ 



Integration with respect to the azimuthal angle gives 

oo 



(r + xY + 2' 



4rx 



(r-x) 



< 0. 



(5) 



As an aside, we remark that the velocity field ii^(r, z) could be 
equally well calculated without the intermediate step of determining 
air). In this respect, we should first find the kernel function 

CX) 

K{r, z-x)= \ AJi (Ax) Ji (Ar) exp (-A \z\) dA, x, > 0. 



In particular, in the limit jz — )■ 0, K(r, z; x) can be expressed by the 

Dirac 5-function, K{r, 0; x) ~ r'^^Sir — x) = (^x)"^''^ ^{^ ~ x)- 
With the aid of this kernel, we can show, on combining equations (Al) 
and (A5) from the appendix in Bratek, Jalocha & Kutschera ( 2008) 
and using equation that the following concise relation between 
the velocity field v^(r,z) and the rotation curve v^[r,Q) = Vcr{r) 
holds: 



{r,z) = r / K{r,z;x)Va{x)dX- 



(6) 



Obviously, condition {4}, and thus equations Q and {6}, hold for az- 
imuthal velocity only approximately, as orbits of test bodies form very 
complicated spatial curves, occasionally crossing the galactic mid- 
plane. These equations, however, should hold approximately for or- 
bits whose projections on to the galactic mid-plane are almost cir- 
cles. For such orbits, equations (O or (|6} could be used for estimating 
the vertical gradients of azimuthal velocity by direct calculation of 
dzV^(r, z). Before proceeding further, we first try to convince our- 
selves that the above approximation can indeed be applied. 

By taking into account the first integral of motion existing by 
axial symmetry, r'^{t)ip'{t) — jz, where jz s a constant characterizing 
a given trajectory, the radial component of the equations of motion 
r"{t)-r{t){(p'{t)f = -9r$(r(t),2(t))canberecastintheformof 
an equation for the image of the trajectory projected on to the galactic 
mid-plane: 



r^{ip) + 2r'{ipf - r{ip)r"{if) _ r{ifi)gr{r{ifi), z{ip)) 



(7) 



Here, v^{r{(p), z{ip)) = r{(p)(p' {t{(p)) = jz/riip) and 
gr{r{Lp), z{(fi)) = —dr(y,)^{r{(f),z{Lp)) is the gravitational accel- 
eration experienced on the trajectory. The term on the left-hand side 
of equation l|7j is proportional to the curvature Kxy{f) of the projec- 
tion 



(r2(^) + r'2((p))=*/' 

Thus, equation ^ is exactly satisfied for orbits for which the left- 
hand side of equation ^ identically equals 1. This is true for a 
trajectory whose projection on to the galactic plane forms a circle, 
r{(p) = const., or a hyperbolic spiral, r{ip) = a{tp — tpo)^^- As for 
realistic trajectories, equation (5) should work in the approximation 




Figure 1. Distribution functions for the radial force equilibrium test (solid 
line) and for the orbit circulaiity test (dotted line); these smooth-looking lines 
were obtained by filtering data using averaging in a window of small x-width. 
The dotted line shows the superposition of distribution functions of variable 
X = logj^Q J measuring departure of radial distance r from its average value 
f on a given orbit, obtained for various trajectories (here, x = 0.03 and 
5x = 0.3). The solid line shows the analogously obtained distribution func- 
tion for variable x = logj^Qt;^/ {rgr), (here, x = —0.014 5x = 0.09) 
measuring the depaiture of centrifugal force' v'^/r from the actual magnitude 
of the radial component of the gravitational acceleration. 



of quasi-circular orbits. A quasi-circular orbit by definition can be en- 
closed within a relatively thin toroidal tube symmetric with respect to 
the galactic mid-plane and concentric with the centre. The projection 
of such orbits on to the galactic mid-plane is roughly circular. 

The expectation that equation {5} holds approximately for quasi- 
circular orbits has been confirmed by the following numerical experi- 
ment. Starting from various initial conditions, trajectories of test bod- 
ies were found by numerical integration of equations of motion in 
the gravitational field of the disc of surface density ([3} correspond- 
ing to the SHO rotation curve. The initial positions were generated 
randomly in space with probability density proportional to the sur- 
face mass density and falling off exponentially with z. At the initial 
instant of each simulation, the only non-zero velocity component as- 
sumed non-zero was the azimuthal one. For a given initial position, 
the initial velocity was calculated from equation (|5]l. To test our hy- 
pothesis, we found nearly 300 trajectories with various complicated 
shapes. For the majority of these, the radial variable on a given or- 
bit was observed to have small dispersion relative to the mean radial 
variable. The quasi-circular orbit condition (|4j was also satisfied with 
small dispersion. The relevant statistical analysis is presented and ex- 
plained in more detail in Fig. [T] From the point of view of the pre- 
cision of astronomical data and the accuracy of rotation curve mod- 
elling, such a result can be regarded at least as satisfactory. This way, 
the approximate relation (O has been justified. 

For our needs, approximation Q is also justified by the fact 
that the results for gradient estimates in iLevine et al.l | |2008|) , which 
we aim to reconstruct in the global disc model approximation, were 
obtained under the idealized assumption of non-intersecting gas or- 
bits with vanishing velocity in the ^-direction (including orbits off the 
galactic mid-plane). This is also the central assumption in the method 
of deriving the azimuthal velocity profile as a function of the radial 
distance from the Galactic Cent r e, at least for the part of rotation 
curve of interest in iLevine et al.l ( l2008l ). Our approximation is thus 
sufficient, as it is no less accurate than the method of obtaining the 
rotation curve on different heights off the Galactic mid-plane. Hence, 
there is no reason or need to make our equations more realistic than 
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Figure 2. Vertical gradient of azimutlial velocity for quasi-circular orbits cal- 
culated in a disc model (thick line), in a three-component model with a light 
dark halo (dashed line), in a three-component model with a massive dark halo 
(dotted line) and in a maximal halo model (thin Hne), all discussed in the 
text, shown for different values of radial variable r g (3, 8) kpc in steps of 
Ar = 0.5 kpc. 



they are in the present form. All in all, what we indeed assume in our 
paper, which may be considered non-standard, is that gross dynamical 
Galaxy mass forms a flattened rather than spheroidal object. 



2.1 Analytical estimates of vertical gradients in the global thin 
disc model 

The vertical gradient of azimuthal velocity for quasi-circular orbits 
can be found directly from equation ([5} by differentiation with respect 
to z: 



dzV^ {r, z) = 



Gz 



(8) 
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Various gradient lines are depicted in Fig. |2] (including the analogous 
lines of other models studied in the following sections). It is seen from 
this figure that the absolute gradient value decreases with the growing 
mass of the spherical mass component. It is also evident that the disc 
model is naturally suited to describe large vertical gradients of rota- 
tion observed close to the galactic mid-plane. A significant spherical 
component reduces the gradient value in two ways: it has a neces- 
sarily small contribution to the resultant gradient close to the galac- 
tic mid-plane; it removes mass from the disc component, reducing in 
turn the disc contribution to the overall gradient. The fact that the disc 
component gives a larger contribution than the spherical component 
to the vertical gradient in the vicinity of the galactic mid-plane is well 
illustrated in the example of the Kuzmin disc discussed in Section|5] 
In what follows, we determine a global vertical gradient in the 
Galaxy by averaging the behaviour of local rotation fall-offs in the 
region r £ (0, 8) kpc in a way that mimicks gradient measure- 
ments in galaxies. This enables a comparison with esti mates of verti- 
cal derivatives of the Galaxy rotation curve reported bv lLevine et al.l 



Figure 3. Variation of azimuthal velocity on quasi-circular orbits as a func- 
tion of distance off the Galactic mid-plane in the global disc model. Each 
point represents the mean azimuthal velocity component over the radial inter- 
val r g (3, 8) kpc at a given height off the Galactic mid-plane, calculated 
using equation (5). The dashed hnes represent the corresponding one-standard 
deviation band. The linear regression fit slope is — 21.3±4.2 km kpc~^. 



( l2008l) . To this end, we first prepared an array of azimuthal veloci- 
ties calculated from equation (O for pairs (r, z), r £ (3, 8) kpc and 
z e (0.22, 2.62) kpc, in steps of Ar = 1 kpc and Az = 0.2 kpc. 
This assumes that the weighting function is homogeneous in space. 
Next, for a given z, we calculated the mean azimuthal velocity in 
the radial interval and the corresponding standard deviation from the 
mean in this interval. The mean velocity is a reasonable estimate, as 
the Galaxy rotation curve is roughly flat over this interval. The global 
vertical gradient and its error can be now determined by finding the 
slope of a linear regression fit to these data (see Fig.[3]l. 

In this way, we find that the vertical gradient of azimuthal veloc- 
ity for quasi-circular orbits in the global disc model approximation is 
—21.3 ± 4.2 km kpc~^, whereas th e Galaxy gradient measure- 
ments yield -22 ± 6 jLevine et al .1120081) . 

The preceding analysis has assumed that the variable z interval is 
greater t han and lies o utside the measurement interval \z\ < 0.1 kpc 
used in jLevine et alBoOS) . However, because of the linearity of the 
fall-off of the azimuthal velocity evident in Fig. [S] the disc model 
gradient value can be extrapolated in the direct vicinity of the galac- 
tic mid-plane, where the actual gradient should be comparable (the 
disc model might give absolute values of the gradient at the galac- 
tic mid-plane greater than those observed, but this would only be 
the model effect caused by the assumed infinitely thin layer of disc 
mass). For more realistic finite width discs, the absolute gradient 
value would be lower at the Galactic mid-plane, but it would still re- 
main slightly greater than or comparable with the gradient values at 
larger heights off the plane, larger than the stellar disc width (of the 
order of 0.3kpc), where both the thin layer disc and a finite width disc 
must surely give comparable predictions. As seen in Fig. [3] the thin 
disc gradients are almost constant at heights several times greater than 
the width of the Galaxy stellar disc. We stress the observational fact 
that the absolute gradient values off the Galactic mid-plane, which 
we predict for the Galaxy, are consistent with the fall-off measured in 
haloes of ot her galaxies (compare the discussion in the summary of 
jLevine eTlil. 2008)), This suggests that the gradients are indeed con- 
stant at large distances off Galactic mid-planes. Because for remote 
galaxies such measurements must necessarily cover a broader z inter- 
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Figure 4. Azimuthal velocity curves, v^{r,z), obtained in the global disc 
model approximation for quasi-circular orbits as a function of distance from 
the Galactic mid-plane, shown in steps of Az = 0.6kpc. The absolute value 
of the vertical gradient of azimuthal velocity grows with the separation Av^ 
of the curves in this diagram. 



val, as we have assumed, this gives observational support for the fact 
that the gradient should be almost constant close to the Galactic disc. 

For brevity, we call the above method of global gradient deter- 
mination the 'I-method'. The I-method works well when calculations 
are carried out over the flat part of the rotation curve, otherwise the 
velocity dispersion might be large, leading to large uncertainties in the 
gradient determination. However, the global gradient in the region of 
interest can also be estimated by finding a mean value over a radial 
interval of local vertical gradients determined at each r separately. 
By a local gradient at a given r we mean a slope of a linear regres- 
sion fit to the values of azimuthal velocities calculated from equation 
l[5]l, by assuming various z in the same region as in the I-method. We 
call this method of global gradient determination the 'Il-method', to 
distinguish it from the I-method. A global gradient obtained thus is 
— 21.08 ± 5.75 km kpc~^. The standard deviation of this value, 
relative to the mean value, measures the degree of change of the gra- 
dient with the radial variable. Fig.|4]illustrates this change by showing 
rotation curves at different heights off the Galactic mid-plane. 

Finally, it is worth noting that the vertical fall-off of the rotation 
curve determined by a linear fit 



r{r) + jz + S\ 







to the data array considered in this section is 7 = 21.2 ± 
2.4 km s^^ kpc^^ with 95 per cent confidence level, again consis- 
tent with the measurements. We have omitted the rolling term as 
the disc model assu mes z-reflection sym metry. This method is analo- 
gous to that used in jLevine et al1l2008h for the determination of the 
gradient value from Galaxy rotation measurements. 



2.1.1 Vertical gradient of rotation close to the Galactic mid-plane 
(at \z\ < 0.1 kpcj 

Above, we have estimated the vertical gradient for the Galaxy using 
the rotation speed (equation l|5]l) outside the strip \z\ < 0.1 kpc. This 
choice was dictated by the desire to avoid difficulties in numerical 
integration close to the mid-plane where integrands in equation l[8) 
become divergent at 2: = and thus numerically intractable close 
to 2 = 0. We have also given some arguments for the fact that our 
gradient estimates could be extrapolated toward the close vicinity of 



Figure 5. Behaviour of the vertical gradient in Galaxy rotation speed at small 
heights off the Galactic mid-plane. The thick line shows the mean gradient 
value at a given 2 averaged over interval r £ (3, 8) kpc, calculated with the 
aid of equation (s). The thin lines show the one-standard deviation from the 
mean value. 



the mid-plane \z\ < 0.1 kpc and thus could also be compared with 
the value determined from observations in this region. However, to 
substantiate our statements, we independently attempted to calculate 
the gradient in the region \z\ < 0.1 kpc directly from integral l[8), 
by carrying out regularized integration in the principal value sense. 
Numerically, this could be achieved by applying the integration rule 



.f{r,x) dx 



f{r,x)dx + 



f{r,x) dx, 



with e being some tiny number, in practice of the order of 10^ . As 
seen from Fig. [5] on the lower height scales \z\ < 0.1 kpc the pre- 
dicted gradient value does not change noticeably. It smoothly overlaps 
with the upper region value and is still consistent with the observed 
gradient value. This is also the upper bound on the absolute gradi- 
ent value expected from other mass models; in particular, a finite disc 
absolute gradient value should be comparable, provided the disc com- 
prises gross dynamical mass. 



2.2 Vertical gradient of rotation from motions of test bodies 

The most convincing estimation of the vertical gradient in the frame- 
work of the global disc model seems to be the statistical analysis of 
trajectories of test bodies. The trajectories are found by numerical in- 
tegration of the equations of motion in the Galaxy gravitational field 
produced by the disc. The vertical velocity gradients are estimated 
by the slope of a linear regression fit on the plane of values {z, v^), 
typical of each trajectory separately. 

Given a trajectory, we can determine its different characteris- 
tics such as the average radial distance, average azimuthal velocity, 
average absolute distance off the Galactic mid-plane, the z-variable 
dispersion, etc. In our context, it seems most appropriate to use tem- 
poral averaging rather than any other averaging with respect to the 
azimuthal angle. For nearly circular orbits in the Galactic plane. 



Let ip{t) describe the time dependence of the angular position on a tra- 
jectory, Aip = ip{T) — v'(O). and let T be the averaging time. A tem- 
poral average of a function u{ip{t)) along the trajectory is defined as 
{u)t = u{ip{t))dt and, similarly, an angular average is («)aip = 
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these averages would be almost the same; however, the results would 
differ from each other for more complicated motions. The problem of 
choosing appropriate averaging methods is general, and most difficult 
to solve in data analysis. Statistical analysis as such is clear, but the 
real problem is to define the appropriate space of events and the prob- 
ability density defined on it. In this respect, symmetry arguments or 
physics hidden behind a particular problem may help. 

By z-reflection symmetry, it is expected that variable z on av- 
erage, should be zero along a given trajectory. The mean absolute 
distance from the Galactic mid-plane, or standard deviation of vari- 
able 2, can be regarded as the simplest measures of a typical distance 
from the mid-plane on a given trajectory. Unfortunately, the notion 
of the typical distance is a matter of convention. For example, for 
a 2-symmetric Gaussian density function, the two typical distances 
are comparable; however, the standard deviation of variable z from 
2 = is about 1.25 times greater than the mean absolute distance 
from 2 = 0. Another problem is the choice of a weighting function. 
The assumption of the temporal averaging method is tantamount to 
the statement that the weight is uniform in time, and each instant of 
time is treated on the same footing. In this case, more probable are 
events that, on average, last longer. For a stationary system (which 
has energy as a constant of motion), the assumption of temporal av- 
eraging is thus justified. However, in the disc approximation also, an- 
gular averaging is justified because of the assumed axial symmetry. 
In the latter case, the weighting function should be taken as homo- 
geneous in the angular variable. Thus, symmetry arguments cannot 
differentiate between two possible weighting functions, at least for 
the totality of all trajectories satisfying various initial conditions dis- 
tributed axially symmetric. Fortunately, the averaging method and the 
choice of distance measure do not drastically influence the predicted 
vertical gradient values, as far as the accuracy of galaxy modelling is 
concerned. 

The results of our analysis are presented and explained in Fig. 
|6] By assuming the mean absolute distance from the Galactic mid- 
plane as the typical distance on a trajectory, we obtain for the gra- 
dient values the following estimates, —24.0 ± 2.1 kms^^ kpc~^ 
and— 17.8±1.1 kms~^ kpc~^ with 95 per cent confidence, respec- 
tively, for temporal averaging and angular averaging along particular 
trajectories. To carry out the estimates, we took into account about 
600 trajectories with equally spaced initial positions and about 300 
additional trajectories already used in Section[2]with the initial posi- 
tions chosen at random with probability proportional to the disc mass 
density. As we have already noted, the gradient estimates will change 
depending on what typical distance off the mid-plane is chosen. For 
example, the absolute values of the gradient should be expected to be 
lower for the maximum distance as the typical distance. Undoubtedly, 
the most reasonable measure of the typical distance is simply the dis- 
persion of the 2 variable on a trajectory. Then, linear fits to temporally 
averaged data on each trajectory give —22.1 ± 1.9 kms~^ kpc~^ 
with 95 per cent confidence level, whereas the angular-averaging give 
— 16.3 ± 1.0 km kpc~^ for an analogous fit and with the same 
confidence level. As expected, the various gradient values are scaling 
with the assumed typical distance. For example, with temporal aver- 
aging, the ratio of {\z\)t to ^((2 - {z)Tf)T is 0.9146 ± 0.0007, 
which is comparable with 0.92 - the inverse of the respective gradient 
values. 

Concluding Section |2l we note that the various, independent 
methods for the vertical gradient estimates presented in this sec- 
tion give results consistent with each other. Also, they all agree 

J'^ u{ip)dLp. On changing variables, the two averages can be com- 
pared with each other. Noting that r(ip{t))ip' (t) = v^(t), it follows that 
Ma^ = ^ /y "(yW) ^"y[fj) dt (u>T.Thus. (u)a^ 7^ (m)t. unless 
1/3' (t) = const. 
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Figure 6. Results of numerical simulations obtained for trajectories of test 
bodies moving in the gravitational field found in the disc model approxima- 
tion based on the SHO Galaxy rotation curve. The top figure shows the mean 
azimuthal velocity (solid circles) as a function of the mean radial distance ob- 
tained for various trajectories (standard deviations from these quantities are 
represented by vertical and horizontal segments), shown together with the real 
rotation curve measurements (empty circles) and with the SHO model rotation 
curve (solid curve). In the central and bottom figures, solid circles represent 
the mean values of the azimuthal velocity, and the respective standard devia- 
tions are represented by vertical segments; both quantities were calculated for 
each trajectory separately. These averages ai'e shown as a function of the mean 
absolute distance from the Galactic mid-plane (central figure) or as a func- 
tion of the standard deviation of the z-variable (the mean z-variable approxi- 
mately equals zero; bottom figure). To obtain these data, temporal weighting 
was used. The solid lines are the linear regression fits to these data, with slopes 
— 24.0±2.1 and —22. lit 1.9 (with 95 per cent confidence level) in the central 
and bottom figures, respectively. For comparison with the slopes, the empty 
circles represent the data copied from Fig. [3] 
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with the observed global fall-off of rotation from the Galactic 
mid-plane of —22 ± 6 kms~^ kpc~^ determined by iLevine et alj 
12008). In particular, this value best agrees with the slope —22.1 ± 
1.9 km kpc~^ of our linear regression fit to a sequence of pairs 
{ \J {{z — {z)t)^)t, {vip)t} for the sequence of the examined trajec- 
tories of test bodies. 



The azimuthal component of velocity on quasi-circular orbits in this 
field is estimated in the same way as in Section |2] from the ap- 
proximated formula {r,z) = —rgr{r,z) which has there been 
shown to be satisfied in a statistical sense. Again, we come to the 
conclusion that the azimuthal velocity can be decomposed such that 
(r, z) = v'l (r, z) + (r, z) + v\ (r, z), where 



3 COMPARISON WITH VERTICAL GRADIENT VALUES 
IN OTHER MODELS OF THE GALAXY 

It is interesting to see how the results of Section |2] compare with 
those of other, customary simple Galaxy models. For this purpose, we 
first consider a three-component Galaxy model consisting of a central 
spherically symmetric bulge, an exponential thin disc and a spherical 
dark halo. Next, we consider a (spherically symmetric) maximal halo 
model. 



3.1 Vertical gradients in a three-component Galaxy model 

The central bulge volume density pb(r) is obtained in a standard way, 
by deprojection of the corresponding surface (or column) mass den- 
sity (Jb{r) with the help of the following Abel relation: 



pb{r) 



1 f Obim 



The surface density is assumed proportional to the empirical bright- 
ness law of de Vaucouleurs ( 1953) which assumes a constant mass- 
to-light ratio in the bulge component: 

The corresponding circular speed in the disc plane can be written in 
analytical form[f| 



vob (r) 



X 



IG-K^K'^ ^ V 40962 / ' Ro 
The galactic disc is assumed to be exponential 



crd(r) = eexp(-2x) , X 



2eo' 



Finding of the velocity on circular orbits in the disc plane is a textbook 
problem with solution 

«od(r) = ^A-KGe^Qx^ (/o(x)^o(x) - A(x)^i(x)), 



2eo' 



The third component is a spherically symmetric (dark) halo with vol- 
ume density 

I 2 

P'^^'^'^- bi + r^ + z^ - 

The rotational velocity on circular orbits is also a textbook result: 



2 = 0. 



VQhir-) = \^ ATvGaobo ~ ~ arctan ^ , 

Hence, the model rotation curve in the galactic mid-plane is Vc{r) — 

\J ^ob + ^od + '^oh- Outside this plane, for z 7^ 0, the radial com- 
ponent gr{r,z) of the gravitational acceleration can be still found. 



A4{x) is expressed by Meijer G-function iGradshtevn & Rvzhikll200oh 



1133 1 53_9 
8'4'S'8'2'8'4' 8 



M{x) = G='^ X 



Vb{r, z) 



RoX 



'"Ob (X) : 



K ^/r^+z^ 
Ro 



Vh{r, z) = 



^Voh 



The azimuthal velocity for the disc component Vd{r, z) has already 
been given in equation in which we should put (j{r) — Od{r) = 
exp {—r/eo). Next, we calculate the vertical gradient of azimuthal 
velocity: 



VbdzVb + VddzVd + VhdzVh 
y/vl +vj + vl 



(9) 



The calculation of gradients for Vb and Vh is easy. For Vd it has already 
been given in equation ([8), where again we put a{r) = ad{r). 

In the case of the three-component model, we examined two 
example fits that accounted for the Galaxy rotation with various 
dark halo masses. The masses of the bulge, disc and dark halo 
are 1.79 x 1O^°M0, 5.90 x lO^^M© and 8.21 x 10^°Mq (to- 



tal 1.59 X IO^Mq) for the model with smaller halo,( k 



e = 



7.67, Ro = 0.606 kpc, Ho = 2.16 x IO^Mqpc" 
7.45 X IO^Mqpc^^ eo = 3.6kpc, ao = 4.1 x 
lO"^M0pc"^, bo = 18 kpc, Rmax = 20 kpc) and 
1.70 X IO^^Mq, 2.41 X 10"^Mq and 1.28 x 10"Mq (to- 
tal 1.69 X IO^Mq) for the model with larger halo.( k = 
7.67, Ro = 0.558 kpc, Ho = 2.41 x 10=* Mqpc'^ 9 = 
5.62 X 10^MQpc~^ eo = 2.62 kpc, ao = 60.6 x 
10"^Mqpc"^, 60 = 3.31 kpc, = 20 kpc) The rotation 

curves of the two models are shown in Fig. [7] The behaviour of the 
vertical gradient (equation ^) of the azimuthal velocity as a function 
of the distance off the Galactic mid-plane is shown for various r in 
Fig.lH 

To estimate the global gradient, we proceed in the same way as 
in Section im for the same averaging region. The averaging I-method 



gives -18.2 ± 3.0 km s"^ kpc"^ and -13.7 ± 3.0 km kpc 
respectively, for models with smaller and larger dark haloes (see Fig. 
[8]l, while the averaging Il-method gives —18.5 ± 5.9 km kpc~^ 
for the smaller and larger haloes, respectively. Fig. |9] illustrates the 
dependence of the vertical gradient on the radial distance in the quasi- 
circular orbit approximation by showing rotation speeds at different 
heights off the Galactic mid-plane. 



3.2 Vertical gradients in the maximal halo model of the Galaxy 

Finally, we examine how a maximum halo model deals with 
explaining the vertical gradients of rotation in the Galaxy. 
The mass function is assumed to be Keplerian, M{r, z) = 

speed. 



where Vn is the observed circular 



By analogy with the disc model symmetry, we assume the quasi- 
circular orbit approximation in the vicinity of the galactic mid-plane. 
The azimuthal component of velocity is then derived from the balance 
condition of the radial component of gravitational and inertial force 
in cylindrical coordinates, r" 
that is. 



_ r GM{R) _ r RvjiR) 



v^{r,z) = —Vc (R) . 
The vertical gradient reads 



R= Vr2+^, 
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Figure 7. Two least-squares fits of a three-component Galaxy model consid- 
ered in the text fitted to the SHO rotation curve (solid lines), with the smaller 
dark halo (top figure) and with the larger dark halo (bottom figure), shown 
together with the unified rotation velocity measurements collected by SHO 
(empty circles). The model rotation curves are decomposed to the central bulge 
(thin line), the exponential disc (dashed line) and to the dark halo (dotted line). 

v.«, = ^ (^v'^iR) - ^) , R = y;^. 

The gradient at different radii is shown in Fig. |2] Again, we deter- 
mine the global gradients as previously in Section |3] The I-method 
gives —8.9 ± 3.8 kms^^ kpc^^, see figure [Tol and the Il-method 
gives —9.7 ± 4.0 km s^^ kpc^^. FigurefTTIillustrates how the verti- 
cal gradient changes with the radial variable. The total Galaxy mass in 
the maximal halo model is 1.7 x IO^^Mq which is 1.5 times greater 
than the total Galaxy mass predicted by the global disc model. 



4 GALAXY NGC 891 

In studies of the vertical gradients of rotation, the case of the Milky 
Way is exceptional, as the measurements have been carried out very 
close to the galactic mid-plane. We have already given some argu- 
ments to suggest that our gradient estimates above the galactic mid- 



Figure 8. The behaviour of the azimuthal velocity on quasi-circular orbits 
as a function of the distance off the Galactic mid-plane, obtained in a three- 
component Galaxy model considered in the text. The top figure shows the 
model with the smaller dark halo, whereas the bottom figure shows the model 
with the larger dark halo. Each point represents an average over the interval 
r € (3, 8) kpc. The dashed lines represent the corresponding one-standard 
deviation band. 

plane could be extrapolated in the direction toward small z. However, 
to have additional support for this statement, we decided to test our 
approach for the edge-on galaxy NGC 891, with one of the deepest 
ever performed H i observations extending out to 22kpc from the 
galaxy disc. Rotational velocity in the halo was observed to decrease 
in the direction normal to the galacti c mid-plane, with a vertical fall- 
off ra te of about —15 km kpc~^( [Oosterloo. Fratemali & SancisH 
l2007h . In agreement wi th the results are observations of diffuse ion - 
ized gas performed by iHeald. Rand. Benjamin & Bershadv 
giving a vertical gradient in azimuthal velocity in the north-east quad- 
rant of the ionized gas halo of about 17.5 ± 5.9kms~^ kpc~^ 
(z € (1.2, 4.8) kpc, r G (4.02, 7.03) kpc), but with no vertical gra- 
dient detected in the south-east quadrant. 

In order to determine the global vertical gradient, we apply an- 
alytical estimates in the framework of the global thin disc model 
of the galaxy in the quasi-circular orbit approximation, as in Sec- 
tion [2] This time, however, we do not study trajectories of test bod- 
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Figure 9. Azimuthal velocity for quasi-circular orbits in a customary thi'ee- 
component model with the smaller dark halo (top) and with the larger dark 
halo (bottom) as a function of the distance from the Galactic mid-plane, shown 
in steps of Az = 0.6kpc. The empty circles represent the measurements of 
Galaxy rotation. 



ies, as we have already seen for the Galaxy that analytical estimates 
are sufficient, giving similar results. We make calculations in the re- 
gion z G (1.2, 4.8) kpc, r £ (4.02, 7.03) kpc. The vertical gra- 
dients obtained for NGC 891 are -19.9 ± 3.0 kms'^ kpc'^ and 
-19.5 ± 1.7kms~^ kpc"^ using the I-method and Il-method, re- 
spectively. These predictions overlap very well within errors with the 
observed vertical gradient values. 

This result gives additional strong support that the thin disc ide- 
alization works well as a tool for studying vertical gradients of az- 
imuthal velocity in flattened galaxies. 



5 SUMMARY AND CONCLUSIONS 

Throughout this paper, we have given several arguments to suggest 
that the global thin disc model of flattened galaxies is naturally suited 
for describing the large vertical gradients of rotation speed observed 
in the neighbourhood of the galactic mid-plane. The gradients were 
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Figure 10. Vaiiation of the azimuthal velocity on quasi-circulai' orbits in func- 
tion of the distance from the galactic midplane in the maximal halo model. 
Each point represents an average over the interval r € (3,8) kpc, the dash 
lines represent the corresponding 1 -standard deviation band. 
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Figure 11. Azimuthal velocity for quasi-circular orbits in function of the dis- 
tance from the galaxy midplane in the maximal halo model shown in steps of 
Az = 0.6kpc. 



estimated in the quasi-circular orbit approximation, which was estab- 
lished in Section|2]to give trustworthy results. 

The gradient values predicted for the Galaxy in this approxi- 
mation agree very well with measurements when the disc comprises 
gross dynamical mass. Independent estimates of the gradient in this 
model (both analytical and from analysing the motion of test bodies) 
give consistent results. The vertical fall-off in the rotational velocity 
in this model is not very dependent on the height, at least out to 3 kpc 
above the mid-plane. This result is consistent with observations of the 
rotation speed in other galaxies. In other words, the distance from the 
mid-plane is not crucial for the gradient determination. Thus, more 
realistic mass models of the Galaxy, such as a finite width disc com- 
prising the whole dynamical mass, should give similar results. We 
have also found that Galaxy models with significant spheroidal com- 
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Figure 12. Azimuthal velocity for quasi-circular orbits in function of the 
distance from the galactic midplane in the global disc model of galaxy 
NGC 891, shown in steps of Az = 0.6kpc (top figure) and vari- 
ation of the azimuthal velocity component with the distance from the 
galaxy midplane (bottom figure). Each point represents an average over 
the interval r G (4, 7)kpc, the dash lines represent the correspond- 
ing 1 -standard deviation band. Rotation curve of this galaxy comes from 
http : //www. ioa . s .u-tokyo . ac . jp/~ sofue) 



ponent are not consistent with the gradient measurements, and the 
discrepancy grows with the mass of the component. 

With the aid of our model, we have also studied the vertical gra- 
dient in NGC 891. The gradients in this galaxy were measured at 
relatively large heights from the galactic disc compared to the mea- 
surements in our Galaxy. Nevertheless, we again obtained predictions 
in accordance with observations. 

Based on all these results, we can hypothesize, contrary to what 
is impUed by dark halo models, that gross mass distribution in our 
Galaxy is more flattened, disc-like, rather than spheroidal. 
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